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Abstract 

We investigate the impact of many-body interactions on group-velocity slowdown achieved via 
Electromagnetically Induced Transparency (EIT) in quantum dots using three different coupling- 
probe schemes (Ladder, V and A , respectively). We find that for all schemes many-body interac- 
tions have an important impact on the slow light properties. In the case of the A and V schemes, 
the minimum required coupling power to achieve slow light is significantly reduced by many-body 
interactions. V type schemes are found to be generally preferable, due to a favorable redistribution 
of carriers in energy space. 
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I. INTRODUCTION 



Quantum dot (QD) based materials are promising for applications exploiting quantum 
coherence phenomena due to their atom-like properties and long dephasing times.- They 
have been proposed to act as active media in devices for controlling the emission pattern of 
phased array antennas^ or in slow light based all-optical buffers.- 

A particular physical effect that can be utilized for generating slow light is electromag- 
netically induced transparency (EIT). EIT refers to an artificially created spectral region of 
transparency in the middle of an absorption line due to the destructive quantum interference 
arising from two transitions in a three-level system.-*^ By virtue of the Kramers- Kronig rela- 
tions such an absorption reduction is accompanied by a large positive slope of the refractive 
index which translates into a reduced group velocity in vicinity of the resonance. 

Very recently, the first experimental studies of EIT in QD systems have been performed^ 
where a coherent absorption dip in a coupling-probe experiment has been observed for an 
optically thin structure. Semiconductor QD based EIT schemes without real carrier exci- 
tations have been studied using models from atomic physics.-'^'^'^ Such EIT configurations 
involve pumping of intraband transitions whose wavelengths lie in the deep infrared, a regime 
for which high intensity laser operation is very difficult. Carrier-exciting schemes using in- 
terband coupling transitions therefore come into play. Recently, a theoretical description 
with the inclusion of many-body effects for a solid state QD EIT A configuration has been 
reported.— 

Using carrier-exciting schemes one in effect addresses two types of quantum coherence 
phenomena, EIT as well as coherent population oscillatioi>^»i^ (CPO). CPO is a four- wave 
mixing effect based on interference between the coupling and probe fields. One should there- 
for keep in mind that schemes involving carrier excitation generally would contain a mixture 
of the two effects. However, CPO can be ruled out by choosing a setup utilizing orthogonal 
polarization directions for the coupling and probe field thus preventing the possibility of 
interference. 

Concerning the EIT effect; an inherent problem of the carrier-exciting scheme is that the 
carriers excited by the coupling field block the transitions via the Pauli blocking factor, effec- 
tively decreasing the strength of the transitions, making such configurations less attractive 
than those that do not excite carriers. In addition, the excited carriers modify the spectral 
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properties of the system via their mutual Coulomb interaction. Such effects cannot be ac- 
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counted for in the non-interacting model. The work presented in Refs. 
the many-body aspects of a carrier exciting A configuration in a transient regime. The study 
of pulse propagation in a semiconductor slow light medium would generally involve solving 
the coupled Maxwell-Bloch equations. However, under certain circumstances an analysis 
of the steady state properties of the semiconductor Bloch equations alone is adequate. In 
the sense that the linear optical response extracted in this limit is directly linked to the 
propagation characteristics of a wavepacket traveling in an optically thick QD system. Non 
carrier-exciting schemes have only been studied with the inclusion of many-body effects in 
this limit on one occasion,— while studies of carrier-exciting schemes have not, to the best 
of our knowledge, been discussed in the literature. The aim of this paper is to present a 
comparison between different EIT schemes, with and without carrier excitation, that can be 
realized in the same dot structure. We study the EIT generated slow light properties of InAs 
QDs by solving the generalized semiconductor Bloch equations (SBE) in the Hartree-Fock 
approximation. The slowdown capabilities of the Ladder, V and A schemes (see Fig. [T]) ob- 
tained in steady state are compared using two models; the atomic model where interactions 
are disregarded, and the interacting model where many-body effects are taken into account. 



II. THEORETICAL MODEL 



The heterostructure under consideration consists of conical InAs quantum dots (radius of 
9 nm and height 3 nm) residing on a 1.2 nm thick wetting layer (WL), sandwiched between 
two slabs of GaAs. The electronic structure is calculated as the solution to the single-band 
Schrodinger equation for the envelope wavefunction in the effective mass approximation.— 
Using effective electron and hole masses mg = 0.067 rriQ, nih = 0.15 rriQ and a conduc- 
tion/valence band offset CBO = 705 meV, VBO = 363 meV, we find six confined hole 
states (labeled |hO) to |h5)) as well as six confined electron states (labeled |eO) to |e5)), 
all doubly degenerate due to spin. Furthermore, the inherent rotational symmetry ensures 
complete degeneracy of the first and second excited- as well as third and fourth excited state 
for both bands. For each band we also find the onset of a continuous set of delocalized states 
extending into the wetting layer. These WL states are treated as plane waves. The resulting 
energy level structure along with three different EIT-schemes are shown in Fig. [H The dipole 
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FIG. 1: Schematic quantum dot level structure and three EIT configurations. The 
frequency of the intense coupling field is denoted Ucoup while the weak probe field is shown 
as ujprobe- For illustrative purposes the figure has not been drawn to scale. 

selection rules allow for the realization of the three archetypical EIT schemes; Ladder, V 
and A. In the two latter cases, the coupling field excites an interband transition, resulting 
in optical pumping of the dot. We model experiments where a continuous wave coupling 
field is irradiating a homogeneous ensemble of QD's. The steady state system response is 
obtained by applying a weak probe pulse with a gaussian envelope at times later than any 
other timescale of the system relative to the onset of the coupling field, such that transient 
effects may be neglected. In order to effectively utilize the slowdown of the light, as in an 
all-optical buffer, the probe field must propagate within the plane of the active medium. 
Assuming the QD's to lie in the x-y plane, we take the probe polarization along the x-axis 
and the propagation direction along the y-axis. To completely rule out CPO effects we need 
orthogonal polarizations of the coupling and probe fields. For the the V and A schemes 
the coupling field is polarized in the y-direction, but in these situations we let the coupling 
field propagate perpendicular to the QD plane so that we can disregard propagation effects, 
e.g. attenuation, in the coupling field. In the Ladder scheme the coupling field connects 
states |eO) and |e5) and is polarized in the growth direction of the QD's (2;-direction) and 
thus traveling in the QD plane. We disregard propagation effects, as the coupling field is 
effectively connecting two empty states, thus rendering the transition transparent. 

The interband dipole moments connecting |eO) and |h5) as well as |e5) and |hO) are 
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non-zero even in the one-band effective mass description. A detailed strain based 8-band 
k ■ p calculation^^ shows however that these dipole moments are roughly a factor of 20 
larger compared to the one-band result. As our focus is on the influence of the many-body 
interactions on the slow down effects we will therefore assume the k ■ p based results for the 
A^eo/is = 10.24 eA and /Xes/io = 10.14 eA dipole moments. The other relevant dipole moments 
are /ieo/io = 15.55 eA and /ie0e5 = 2.79 eA. 

The linear optical response to the probe, i.e. the susceptibility x{^)^ whose real and 
imaginary part are related to refraction and absorption, respectively, is found from the 
macroscopic polarization P{uj) as x{^) = eoE^luj) i where eo is the vacuum permittivity and 
Ep{uj) is the amplitude of the probe field. 

The time resolved macroscopic polarization component in the direction of the probe field, 
P(t), is computed from the microscopic polarizations according to semiclassical theory:— 

P(t) = i(iVdot$^/i.,P.,(t) + ^$^[/XkPk(t) + c.c.)]). (1) 

id k 

Where Pij and are microscopic polarization components of localized dot states («, j) and 
diagonal interband polarization components of delocalized WL states (k), respectively. In 
this treatment we disregard polarization components relating to transitions connecting dot 
and WL states. Dipole matrix elements between localized states are denoted /ijj, whereas 
/ik is the dipole moment relating to WL states. A^dot is the two-dimensional density of the 
dots in the WL plane, A is the normalization area of the WL, and w is the thickness of the 
active region. 

The microscopic polarizations are the off-diagonal components '^u^u2 i^i 7^ ^2) of the 
reduced density matrix PuiU2^ where z/ refers to either a QD state i or a WL state k. The 
time development of the polarizations are found by solving the SBE in the Hartree-Fock 
approximation, see e.g. Ref.[l8j, given (in the electron-electron picture for the sake of brevity) 
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by 



d 

= ihS^^^^t) ^ -ih'jd'^^^^^, (2) 

d 



where 



^u(t) ~^ ^ ^ {^UU4U3U '^UU4UU-i\ Pu3U4,(t^ (4) 



1/31/4 



are the Hartree-Fock renormahzed single particle energy and generalized Rabi frequency, 
respectively, riy is the diagonal component of the density matrix, i.e. p^y. The term 
—GpuiU2E{t) is the electromagnetic field interaction in the dipole approximation, and the 
matrix elements of the Coulomb interaction are V^^y^y^^^ = f $*^(r)$*2(r')e^/eoeb|r — 
r'\^i,.^{r')^y^{r) d^rd^r' . Coulomb elements are found by approximating the numerically 



evaluated localized dot states 
culated following Refs. 
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w those of a harmonic oscillator; and ^ojy,,^, is then cal- 
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For the situations considered here, screening effects are 
disregarded due to low WL densities. 

Off-diagonal scattering terms S^-^^^it) are approximated by a temperature dependent 
effective dephasing rate 7^; 

Sy,U2{t) ~ --fd'^uiU2- (6) 

Diagonal terms representing collision induced particle exchange processes, are mimicked by 
a nonradiative recombination and a population relaxation towards quasi-equilibrium Fermi- 
Dirac functions fu"^^ , determined by the charge- carrier density and temperature. The scatter- 
ing rates are denoted 7c-c and 7c-p representing carrier-carrier and carrier-phonon scattering, 
and the recombination rate is called 7^^; 

-lc-vVv4 {t) - fui {pi, Ti)] . (7) 



Here fi and T are the chemical potential and temperature of either the plasma (p) or lattice 
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(Z), which are found following the procedure presented in Ref. 
interacting (atomic) model by taking the limit where all Coulomb elements and population 
scattering rates are set to zero. The results presented here use a dot density of N^^ot = 
5 ■ 10^^ m~^, a discretization of the WL into 100 A;-points and a fixed lattice temperature 
of 200 K, for which the literature-'-^ gives scattering rates around 7^ = 1.5 ■ 10^^ s~^, 
7c-c = 2.0 ■ 10^2 g-i^ = 2.0 ■ 10" s'\ and 7„^ = 1.0 ■ 10^ s'^ 



III. RESULTS 



The optical response for the three different schemes using a coupling field with an intensity 
of 26 MW/cm^ is shown in Fig. [2] illustrating both the interacting and the non- interacting 
cases. An immediate difference between the atomic and many-body approach is the change 
in probe field energy towards negative detuning. This is due to the excitonic shift of the 
various probe transitions. An apparent feature of the Ladder scheme (Fig. [2al) is that the 
peaks of the imaginary part of the susceptibility are highest for the many-body model. 
While the distance between the peaks remains the same one can readily see that a larger 
area is covered by the many-body spectrum. This is evidence that oscillator strength has 
been shifted into the |eO) — \hO) transition, in other words it has been Coulomb enhanced. 
Considering the the real part of the susceptibility, the enhancement of the probe transition 
results in more pronounced features, notably a larger slope at zero detuning. Both curves 
are generally shifted upwards from the zero point; this behavior is due to transitions in the 
vicinity of the probe which are contributing to the background index of the area in question. 

The optical responses for V and A schemes are shown in Fig. [2b] and Fig. [2cl respectively. 
For the interacting case the coupling field has been tuned to the zero density exciton res- 
onance of the probe transition. The asymmetry in the peak heights of the imaginary part 
of the susceptibility has different origin for the interacting and noninteracting cases. The 
skewness in the atomic model is due to the fact that we are not dealing with a closed three 
level system. The control field pumping the |eO) — \ hO) is also connecting the dipole allowed 
|e5) — \h5) transition, however severely negatively detuned. Effectively we are dealing with 
two EIT schemes, the original V (A) and a detuned A (V) scheme. In general a negatively 
detuned A or \^ EIT scheme has a prominent shift in peak height towards positive detunings 



7 



0) 



a: 



0) 

a: 



10 

8 
6 
4 
2 

-2 
-4 
-6 
-8 
-10 



1 


1 „ 1 


1 1 
W.o. m-b 


- 




W. m-b 

- 








1 


1 1 


1 1 



1 1 r 

W.o. m-b 

W. m-b 



-10 



1 

0.5 


-0.5 
-1 
-1.5 



1 


1 


1 1 1 1 1 
W.o. m-b 






/ W.ji 


1-b 














/ A 








y y^ 1 




/ / 








1 




r 1 ' 
1 1 1 1 


1 


-25 - 


20 


-15 -10 -5 5 



10 

a) Ladder scheme. 
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FIG. 2: Ladder, V and A schemes using a coupling field intensity of 26 MW/cm^: 

Complex susceptibility vs. probe energy. The dashed line corresponds to the atomic 

model, i.e., without many-body interactions, while the solid line is evaluated with 

many-body interactions. The energy is measured relative to the zero density excitonic 

resonances c^q^'^q, e^g^^Q and e^.Q^'^g, respectively. To facilitate a comparison between the two 
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cases, the noninteracting spectra have been shifted accordingly. 



(resembling an optical Stark shift). What is seen in the two figures is an admixture of the 
symmetric peaks owing to the resonant V (A) scheme and a Stark shifted transition. Further 
evidence of this effect has been obtained by altogether disallowing the "conflicting" transi- 
tion in which case one recovers the symmetric result. In the models including interactions 
the same asymmetry should be expected, but it is countered by a negative shift in resonance 
energy induced by the Coulomb interaction with the excited carriers. This means that the 
coupling is detuned positively with respect to the resonance, and hence the asymmetry tends 
towards negative probe energy. The probe transitions in these cases are not enhanced rather 
they are suppressed by the inclusion of many-body effects; this can be seen by the fact that 
the features are generally smaller in magnitude than in the atomic model. The splitting of 
the peaks is larger though, showing that the effective Rabi frequency is higher, owing to the 
Coulomb enhancement of the |eO) — \hO) coupling transition. 

As a basis for comparing the slow light capabilities of the different schemes we examine 
the maximum obtainable slowdown factor S, which is equal to the group index and is a 
measure of the group velocity reduction. The slowdown factor 5* is a flgure of merit relevant 
for optical storage, and is defined via 

e Co On 

S = — = n + uj— , (8) 

Vg OUJ 

where Cq is the speed of light in vacuum, and n = Re{[n^ + x(co')]^} is the refractive index. 
The maximum slowdown is found at the frequency for which the slope of the refractive 
index is largest. Notice, that the slowdown factor obtained away from resonance is given 
by the background refractive index. To make a just comparison, we detune the coupling 
field used in the many-body V and A models from the zero density exciton resonance, so 
that the peaks of Im{x(ti;)} become symmetric. However, the amount of detuning for this 
to be realized depends on the intensity of the coupling field. Changing the intensity in turn 
changes the amount of carriers being excited and thus the excitonic shift, which results in 
different detunings for different intensities. 

Figure [3] depicts the maximum achievable slowdown factor and the corresponding absorp- 
tion coefficient a{uj) = 2^Im{n(u;)} as a function of coupling power for all three schemes. A 
striking feature of Fig. [3] is that the inclusion of many-body effects leads to different results 
depending on the choice of EIT scheme. The results of the carrier-exciting V and A setups 
are seen to differ fundamentally from the Ladder scheme. Inspecting the absorption coef- 
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FIG. 3: (Color online) Maximum slowdown factor and corresponding absorption coefficient 
vs. coupling intensity for various EIT-schemes. The dashed curves are for the independent 
particle (atomic) model, while the solid curves include many-body interactions. 

ficient plot we see two plateaus, corresponding to the maximum (minimum) absorption in 
the absence (presence) of EIT. The transition from the upper to the lower plateau happens 
across fewer orders of magnitude in the coupling power for the ladder scheme than the other 
two schemes. Here the absorption coefficient drop is solely due to the quantum coherence 
effect setting in, whereas the transition for A and V happens across a significantly larger 
relative range. At low coupling power the absorption drop is driven by the excitation of 
carriers occupying the probe transition. As evidenced by the slowdown plot, the EIT effect 
sets in at larger coupling powers, only when we are near the lower plateau. The largest 
slowdown values are achieved using the Ladder scheme, for which the slowdown factor is 
increased significantly when interactions are included. This is due to Coulomb enhancement 
of the |eO) — \hO) resonance probed in this scheme. The slowdown effect is seen to disap- 
pear at the same value of coupling power for both cases, which indicates that the coupling 
transition (|eO) — |e5)) utilized in the Ladder configuration is unchanged by the inclusion 
of many-body effects. On contrary, for the carrier-exciting schemes V and A, many-body 
effects have a significant impact on the coupling threshold. 

Both V and A schemes show largest slowdown values for the non-interacting model. On 
the other hand, the noninteracting model overestimates the minimum required coupling 
power for observing slow light by roughly a factor of two compared to the more realistic 
case of interacting particles. As both schemes utilize the same coupling transition they 
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experience the same coupling power threshold, in both the interacting as well as the non- 
interacting case. The shift in required coupling power can be attributed solely to the Rabi 
energy enhancement of the coupling transition. This conclusion is reached by inspecting the 
absorption spectrum in absence of a coupling field. By comparing the height of the |eO) — \hO) 
resonance (coupling transition used in both schemes) with and without interactions we find 
that the dipole moment of the transition is enlarged by roughly a factor of 1.2. If we, in the 
atomic model, enlarge the coupling dipole moment by the same amount, we end up with a 
result having the same minimal requirement on coupling power as the interacting case. This 



result stands in contrast to the findings in Ref. uM, where a shift in required coupling power, 
due to Coulomb enhancement, of two orders of magnitude was reported. However, this work 
was performed in a transient regime and a direct comparison is therefore not applicable here. 

The V scheme is preferable to the A scheme, due to its higher slowdown values. The 
reason is twofold. Firstly, based on observations from the absorption spectrum without 
coupling, we find that the two probe transitions are both Coulomb suppressed, however the 
A scheme to a higher degree than V. Secondly the fact that the V probe connects a hole 
ground state to an electron excited state results in a larger Pauli blocking factor (the third 
term in equation ([2])) as compared to the A scheme. Figure H] demonstrates this for the 
interacting case utilizing a coupling intensity of 2.5 MW/cm^. For the interacting model 
the redistribution of carriers plays a crucial role; the smaller effective electron mass leads 
to a larger energy spacing of the electronic levels, which means that the electron excited 
states become less populated than their hole counterparts. In the case presented here the 
hole ground state is seen to be depleted and carriers are redistributed into the higher lying 
energy states. For electrons the redistribution is less prominent, and as the coupling field 
excites more and more carriers, electrons accumulate in the ground state. Thus the Pauli 
blocking factor seen by the probe in the V configuration is always the larger, which ultimately 
translates into an increased slowdown factor. This result is quite general and could act as a 
pointer for experimental realization of KIT mediated slow light. 

IV. CONCLUSION 

In conclusion, we have investigated the slow light properties of InAs QDs using a 
model including many-body effects for three different KIT schemes and found fundamental 
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FIG. 4: Temporal development of the Pauli probe blocking factors for A and V schemes at 

coupling intensity 2.5 MW/cm^. Shown in the figure are the Rabi oscillations of the 
populations at the onset of the coupling field and the relaxation towards their stationary 
value. Also shown is an illustration of the "turn on" of the coupling laser. 

differences. The Ladder scheme that utilizes a transparent coupling transition has its 
slowdown factor increased due to Coulomb enhancement. However, there is observed no 
change in the necessary coupling power required to reach EIT. Conversely in the V and A 
schemes, many-body effects enhance the coupling transition resulting in a lowering of the 
necessary coupling power. The V type configuration is found to be preferable, due to a 
favorable redistribution of carriers. 
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